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ABSTRACT 

We compute CBR anisotropies in mixed models with different hot components, 
including neutrinos or volatile HDM arising from the decay of heavier particles. The 
CBR power spectra of these models exhibit a higher doppler peak than CDM, and 
the discrepancy is even stronger in volatile models when the decay gives rise also to a 
neutral scalar. 

CBR experiments, together with Large Scale Structure (LSS) data, are then used 
to constrain the space parameter of mixed models, when values of the primeval spectral 
index n > 1 are also considered. Even if n > 1 is allowed, however, LSS alone prescribes 
that 0.30. 

LSS can be fitted by taking simultaneously a low derelativization redshift Zder 
(down to ~ 600) and a high n, while CBR data from baloon-borne experiment cause a 
severe selection on this part of the parameter space. In fact, while late derelativization 
and 71 > 1 have opposite effects on the fluctuation spectrum P{k), they sum their 
action on the angular spectrum C/. Henceforth n> 1.3 seems excluded by baloon- 
borne experiment outputs, while a good fit of almost all CBR and LSS data is found 
for fi/i values between 0.11 and 0.16, n ~ 1.1 and z^er ^2000-5000. A smaller n is 
allowed, but Zder should never be smaller than ~ 1200. 

Key words: dark matterrdecaying particles, dark matter: massive neutrinos, large 
scale structure of the Universe, cosmic microwave background: anisotropies. 



1 INTRODUCTION 

Anisotropies in the cosmic background radiation (CBR) are 
a strong potential source of information on the cosmologi- 
cal model. Unfortunately, anisotropy observations are hard 
and significant measures were obtained only recently. As a 
matter of fact, the theory of CBR anisotropies is well under- 
stood (see, e.g., Hu & Sugiyama 1995 and references therein) 
and public numerical codes allow to calculate the expected 
anisotropies for a wide range of cases (Seljak & Zaldarriaga 
1996). It is then easy to see that CBR anisotropies depend 
on all the ingredients that define a cosmological model: the 
background metric, the substance mix and the primeval fluc- 
tuation spectrum. Several authors used available codes to 
predict CBR features for suitable ranges of model param- 
eters. However, within the range of models consistent with 
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the inflationary paradigm, not enough attention, in our opin- 
ion, has been devoted yet to mixed models. Anisotropies 
expected for them were calculated by Ma & Betschinger 
(1995), De Gasperis et al. (1995), Dodelson et al. (1996), 
but parameter choices were restricted to cases for which 
anisotropies only marginally differ from the standard CDM 
case. Here we plan to extend the analysis to a wider set of 
mixed models including those for which a greater discrep- 
ancy from standard CDM can be expected and, in partic- 
ular, models with primeval spectral index n > 1 and late 
derelativization of the hot component. If hot dark matter 
(HDM) is made of massive i/'s, there is a precise constraint 
between its density parameter (Qh) and its derelativization 
redshift (zder)'- 



^ 0.68(2der/10*)((;,/6) 



(1.1) 



(see eqs. 2.3-2.4 below; g^, is the number of i/ spin states). 
Henceforth, in order to have 0.10-0.15, Zder cannot be 
lower than ^ 2000, even for g^, — 6. In order to have lower 



© 0000 RAS 



2 E. Pierpaoli and S. Bonometto 



Zder and/or greater Q.^, HDM must arise from the decay 
products of heavier particles. In fact, decay products have 
extra kinetic energy arising from mother particle mass en- 
ergy and therefore have a later Zder ■ Several authors consid- 
ered such scenario, assuming the decay of a heavier neutrino 
into a lighter one (Bond & Efstathiou 1991, Dodelson et al. 
1994, White et al. 1995, McNally & Peacock 1995, Ghizzardi 
& Bonometto 1996), and there are recent attempts of con- 
straining it using CBR data (Hannestad 1998). However, a 
wider range of models give rise to similar pictures, e.g., if 
metastable supersymmetric particles decay into lighter ones 
(Bonometto et al., 1994; Borgani et al. 1996). In a num- 
ber of recent papers HDM arising from decays was called 
volatile, to stress its capacity to grant a later derelativiza- 
tion, which weakens its contribution to the formation of in- 
homogeneities. 

CBR anisotropics were first detected by the COBE- 
DMR experiment (Smoot et al. 1992). The angular scales 
observed by COBE were rather wide (~ 7°) and allowed 
to inspect a part of the spectrum which is almost substance 
independent. Nevertheless, COBE measurement provide the 
normalization of density fluctuations out of the horizon, and 
fair constraints on the primeval spectral index n (Bcnnc>t et 
al. 1996). More recent baloon-borne and ground based ex- 
periments investigated CBR fluctuations on scales compaxa^ 
ble with the horizon scale at recombination In the standard 
CDM model, on these scales, one expects the first doppler 
peak and, unlike what happens on larger scales, anisotropics 
are related both to the spectral index n and to substance 
mix. In principle, these scale are the best ones to test mixed 
models, as on even smaller scales (/> 500 see below) the 
CBR spectrum can be distorted by other effects, like reion- 
ization and lensing. 

The degree scale measurements currently available seem 
to have detected the doppler peak at a fair angular scale, 
but with an amplitude higher than expected in a standard 
CDM scenario with n = 1 (Scott et al. 1995, de Bernardis 
et al. 1997). An analysis of such outputs seems to exclude 
low values of the density parameter f2 = p/ per (Hancock et 
al. 1998, Bartlett et al. 1998). Here p is the present average 
density in the Universe and pcr = SH^/SttG depends on 
the value of the Hubble parameter H = h lOOkms^^Mpc"^ 
Furthermore, Lineweaver (1997) and Lineweaver & Barbosa 
(1998) outline that observed anisotropics are still too large 
to agree with CDM and n = 1 unless h ~ 0.3. Alternative 
causes for such greater fluctuations can be n > 1 or a cosmic 
substance comprising a substantial non-CDM component. 

A number of large scale structure (LSS) observables are 
obtainable from the linear theory of fluctuation growing. A 
wide range of mixed models predict fair values for them. 
Tests of non linear features were also performed through 
N-body simulations (see e.g. Ghigna et al. (1997) and ref- 
erencies therein) . Although some technical aspects of mixed 
model simulations arc still questionable, one can state that 
suitable DM mixtures allow to fit LSS data from 1 to 100 
Mpc, almost up to scales covered by the current CBR exper- 
iments, so that a simultaneous analysis of CBR anisotropics 
and LSS allow complementary tests of the models. 

It ought to be outlined that n > 1 and a hot compo- 
nent have a partially compensating effect on LSS and, in a 
previous paper (Bonometto & Pierpaoli 1998, hereafter BP; 
see also Lucchin et al. 1996, Liddle et al. 1996), we discussed 



how they could be combined to obtain a fit with LSS data. 
On CBR fluctuations, instead, they add their effects and a 
quantitative analysis is needed to see how far one can go 
from both n = 1 and pure CDM. Most of this paper is fo- 
cused on such analysis and on the tools needed to perform 
it. 

In particular, let us outline that public codes, like CMB- 
FAST, cannot predict CBR anisotropics for mixed models 
with a hot component of non-thermal origin. A part of this 
work, therefore, required a suitable improvement of current 
algorithms. 

In this work only models with A = 0, /i = 0.5 and f2 = 1 

are considered. The cosmic substance is fixed by the partial 
density parameters: Oi, = pb/pcr for baryons, fie = pc/ per 
for cold-dark-matter (CDM), Qh and Zder (sec below) for 
HDM. Here we shall also distinguish between HDM made of 
massive f's {neutrino models) and HDM arising from heav- 
ier particle decay {volatile models). The relation between 
the nature of HDM and the amount of sterile massless com- 
ponents (SMLC hereafter) needs also to be discussed (in 
standard CDM, SMLC is made by 3 massless i^'s). 

Early deviations from homogeneity are described by the 
spectrum 

P*(fc) = A^x^{xky-'' (1.1) 

{x = Xo — Xrec IS thc distaucc from the recombination band, 
as Xo is the present horizon radius). Here k — 2n/L and 
L is the comoving length-scale. Models with 1 < n < 1.4 
were considered. In Appendix A we review which kinds of 
inflationary models are consistent with such n interval. 

Section 2 is dedicated to a brief discussion on the differ- 
ent kinds of hot dark matter that could lead to a mixed dark 
matter scenario. In section 3 we analyze the CBR spectrum 
of such models, distinguishing between the effects due to the 
SMLC and those due to the actual phase-space distributions 
of the hot particles and discussing how current algorithms 
need to be modified to provide CBR anisotropics for volatile 
models. We then perform an analysis of the parameter space: 
models are preselected according to LSS constraints related 
to the linear theory. This selection is based on BP results, 
whose criteria will be briefly reviewed (section 4). BP re- 
sults, however, were restricted to the case = 0.1. Here 
we shall inspect a greater portion of the parameter space, 
by considering also models with fi^ = 0.05, allowing for a 
substantial dependence of the CBR power spectrum on the 
baryon abundance. Section 5 is dedicated to a comparison 
of the CBR spectra with current available data, and to the 
final discussion. 



2 DARK MATTER MIX 

The substance of mixed models can be classified according 
to their behaviour when galactic scales (lO^-lO^^M©) enter 
the horizon. Particles already non-relativistic are said cold. 
Their individual masses or energy distributions do not affect 
cosmological observables. Tepid or warm components (see, 
e.g. Pierpaoli et al., 1998) become non-relativistic while 
galactic scales enter the horizon. Hot component(s), instead, 
become non-relativistic after the latter scale has entered the 
horizon. Neutrinos, if massive, axe a typical hot component. 
They were coupled to radiation for T > Tv,dg ~ 900 keV. 
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If their mass m ^ T,^,dg, their number density, at any 
T < T^^dg, is = (3C(3)/47r2)g^T^3 (after electron anni- 
hilation = T^,da[ai.,dg/a{t)] = (4/ll)^/^T, where T is 
radiation temperature) and their momentum distribution 
(normalized to unity) reads 

also when p <^ m. Henceforth, when T <^ m, their distri- 
bution is not thermal, although its shape was originated in 
thermal equilibrium. Notice that, for high p, $1, is cut off as 
exp(-p/T.). 

Using the distribution (2.1) we can evaluate 

(p) = (77rVl80C(3))r, = 3.152 . (2.2) 

If we define Zder as the redshift for which (p) = m, eq. (2.2) 
tells us that Zder occurs when T^, = 0.317 m. In the following 

we shall use the parameter d = lO^/zder, which normalizes 
Zder at a value (iC) in its expected range. At z = 10*, 
photons and CDM would have an equal density in a pure 
CDM model with h = 0.5 and a present CBR temperature 
To = 2.7333 K. Henceforth, such redshift is in the range 
where wc expect that rolativistic and non-rolativistic com- 
ponents have equal density (equivalence redshift; Zeq); be- 
sides of photons and SMLC, it is possible that HDM con- 
tributes to the rolativistic component at Zeq- Its value, in 
different models, is given by cq. (2.11), herebelow. In gen- 
eral, wo shall normalize To at the above value (which is well 
inside 1 a in respect to data) and define 9 = ro/2.7333K. 
For neutrinos of mass m, 

d = 5.2972 • 4h^/{m/eV) , (2.3) 

while 

Qh = 2.1437 • W~^g^{m/eV)9^/4h'^ . (2.4) 

Let us now compare these features with those of a 
volatile model, where the hot component originates in the 
decay of heavier particles X (mass mx), decoupled since 
some early time tdg (see also Pierpaoli & Bonometto 1995, 
and Pierpaoli et al. 1997). If the temperature Tdg -C mx, 
such hot component may have a number density much 
smaller than massive neutrinos. Let Nx,dg be the heavy par- 
ticle comoving number density at decoupling. At t ^ tdg 
their comoving number density reads: 

Nx(t) = Nx,dgcxp[-(t-tdg)/rdy] (2.5) 

with tdg <C Tdy (decay time). Assuming a two-body decay 
process X —>■ v + 4>, into a light (voJatiJe) particle v (mass 
m <^ mx) and a massless particle (j>, it is shown that the 
volatile distribution, at t ^ Tdy, reads 

$„(p,t) = 2(g/p)exp(-g) , (2.6) 

where 

Q = pV/ and p = {mx/2)[ady/a{t)] (2.7) 

provided that X's, before they decay, never attain a den- 
sity exceeding relativistic components causing a temporary 
matter dominated expansion. 

At high p, the distribution (2.6) is cutoff oc 
exp[— (p^/p^)] and this is true also if this temporary regime 
occurs. 

In BP it is shown that, if the massless particle is a 



photon (7), such temporary expansion can never occur in 
physically relevant cases, which must however satisfy the 
restriction il^d <S 1. This limitation does not hold if (f) is 
a massless scalar as is expected to exist in theories where 
a global invarianco is broken below a suitable energy scale 
(examples of such particles are familons and majorons) . 

Using the distribution (2.6), it is easy to see that the 
average 

(p) = iV^/4)mxady/a{t) = p^/2 (2.8) 

and «'s will therefore become non relativistic when p = 
2m/v^; henceforth 

p = 2 • W^zmd/^yn (2.9) 

can be used in the distribution (2.6), instead of eq. (2.7). 

If <j)'s are sterile scalars and the decay takes place after 
BENS, they will contribute to SMLC and affect CBR and 
LSS just as extra massless neutrino states. 

Let us recall that, in the absence of X decay, the ratio 
Pv/ Pt = 0.68132(pi,/6) = Wo- X decay modifies it, turning 
Qv into an effective value 

9^,eff =g. + (16/7)(ll/4)*/3n,,d (2.10) 

In particular, (^'s lower the equivalence redshift. For d(l — 
^h) > 1 + Wo, also v's are still relativistic at equivalence. 
Accordingly, the equivalence occurs at either 

_ 4fe^ 10* _ _ 4/i^ 10*(l-nfe) 

^"•^ ~ 1 + Wo + 0.hd ~ 1 + Wo + 2Q.hd 

(2.11) 

in the former and latter case, respectively. 

Volatile models, as well as neutrino models, can be 
parametrized through the values of Vlh and d. However, at 
given Q,h, the latter ones are allowed only for discrete d val- 
ues (notice that such d values are independent of h and can 
be only marginally shifted by changing 6). The former ones, 
instead, are allowed for a continuous set of d values. This 
can be seen in fig. 1, which is taken from BP, where more 
details on volatile models can be found. In fig. 1 we also 
show which models are consistent with LSS constraints and 
COBE quadrupolo data, for various n > 1. Such constraints 
will be briefly discussed in the next section. In general, they 
are fulfilled for a part of the allowed Q values. Fig. 1 also 
shows that there is a large deal of mixed models with low 
Zder which are allowed by LSS data and are not consistent 
with HDM made of massive i^'s. 

In this article we shall show that a portion of this ex- 
tra parameter space seems however forbidden by CBR con- 
straints. 



3 THE RADIATION POWER SPECTRUM 

To describe the evolution of radiation anisotropics in an ex- 
panding Universe, it is convenient to write the metric in the 
form 

ds^ = a{T)[{l + 2*)dT^ - (1 - 2^)^o,0dx"dx'^] . (3.1) 

in the conformal newtonian gauge. Hero x" (components of 
the vector x) are space coordinates, t is the conformal time, 
a(T) is the scale factor and gives the spatial part of 
the metric tensor in the homogeneity limit. The deviation 
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Figure 1. Parameter constraints from LSS. Dashed lines refer 
to neutrino cases with 1, 2 and 3 massive neutrinos (from left 
to right). Solid lines are the contours of the regions allowed in 
volatile models for different values of n ( marked next to the ar- 
eas). The plot refers to Qj, = 0.1. The constraints considered are 
the shape parameter F, the eg and the number density of clusters 
are considered (see BP for a detailed discussion). 



from a pure Friedmann metric, due to gravitational field 
inhomogeneities, are given by the potentials ^ and 

In the presence of inhomogeneities, the temperature of 
radiation T{n) = {T)[l + Arin)] contains an anisotropy 
term, which can be thought as a superposition of plain waves 
of wave-numbers k. In respect to a given direction h, the am- 
plitude of the single k mode can be expanded into spherical 
harmonics. For our statistical aims it is however sufficient 
to consider the anisotropy as a function oi jj, = k - h and use 
the expansion: 

l — oo 

AT(fc,n,T) = ^H)'(2; + l)A,(fe,r)P,(M) , (3.2) 

1=0 

where Pi are Legendre polynomials and whose coefficients 
can be used to work out the angular fluctuation spectrum 



Ci = {Any / dkk^P^{k)\Ai{k,T)\^ 



(3.3) 



which, for a gaussian random field, completely describes an- 
gular anisotropics. 

At the present time To and for a comoving scale given 
by the wavenumber k, we can compute A; performing a time 
integral (Seljak and Zaldarriaga, 1996) 

Ai{k,To)= [ dxS{k,r)ji[iTo-T)k] (3.4) 
Jo 

over the source function S, which depends upon inhomo- 



geneity evolution inside the last scattering band and from it 
to now. 

The physics of microwave background anisotropies due 
to adiabatic perturbations has been deeply investigated in 
the last few years. 

It has been shown that the characteristics of the peaks 
in the Ci spectrum are related to the physics of acoustic 
oscillations of baryons and radiation between the entry of a 
scale in the horizon and the last scattering band, and on the 
history of photons from last scattering surface to us. 

Background features, like the overall matter and radi- 
ation density content, h and A, have an influence both on 
the positions of the peaks and on their amplitude, but the 
latter also depends greatly on the baryon content Qi, and 
more slightly on the characteristic of the hot component. 

In the following we shall analyze in detail the angular 
spectrum of volatile models, outlining its peculiarities with 
respect to standard CDM and neutrino models. In order to 
do so, we need to modify available public codes, like CMB- 
FAST, allowing them to deal with a hot component whose 
momentum distribution is (2.6). 

It should also be recalled that volatile and neutrino 
models, for given ilh and d, are expected to include a differ- 
ent amount of SMLC. In neutrino models SMLC is less than 
in pure CDM and even vanishes if all v's are massive (unless 
extra SMLC is added ad hoc). In volatile models, instead, 
SMLC is however more than in pure CDM, as scalar <(>'s are 
added on top of standard massless i^'s. 

Several C; spectra of volatile models are presented in 
figs. 6-15. They show two main features, if compared with 
standard CDM: the first doppler peak is higher and the sec- 
ond and third doppler peaks are slightly shifted to the right. 

In principle, we expect volatile model spectra to dif- 
fer from neutrino model spectra because of the momentum 
distribution of volatiles and the extra SMLC they have to 
include. In the following, we shall try to disentangle these 
two effects. 

To this aim we coupled each volatile models with a 
technical neutrino case with identical flh and d, but a 
greater number of neutrino degrees of freedom, so to en- 
sure equal high-redsliift energy densities. In fig. 2 we report 
the scale-factor dependence of the energy densities p{a) of 
volatiles and a massive neutrinos in two coupled models. In 
the case shown, the two energy densities never differ in ratio 
more than 10~^; for different choice of the parameters the 
curve is just shifted to higher or lower redshifts according to 
the value of Zder- 

More in detail fig. 2 states that volatiles have a slower 
derelativization than neutrinos: the transition phase from 
the relativistic to the non-relativistic regime starts earlier 
and goes on for a longer time. This behaviour is related to 
the different shapes of the two distribution functions, and 
to the fact that the volatile one is smoother around (p), 
which corresponds to a value significantly smaller than its 
maximum, after which it is rapidly cutoff [see eqs. (2.6), 
(2.7)]. 

Friedman equations show that T{a)^yp{a} is approxi- 
mately constant. Hence, once we know p(a), we can per- 
form a comparison between the eonformal times of coupled 
volatile and technical neutrino cases. It shows a marginal 
discrepancy as already the p(o) in the volatile and technical 
neutrino cases are very similar, and moreover the hot com- 
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Figure 2. density evolution of a volatile model in comparison 
with the corresponding technical neutrino case. 



Figure 3. comparison of conformal time in a CDM standard 
model and in a volatile model. At high redshifts, volatile models 
typically show a smaller conformal time. 



ponent always contributes as a small fraction of the total 
energy density. On the contrary, if a similar comparison is 
performed between standard CDM and volatile models, big 
discrepancies are found, especially at high redshifts. In fact, 
in the volatile cases the relativistic background is greater 
due to the contribution of the sterile component, and the 
conformal time is therefore smaller than in the CDM case 
(see fig. 3). This implies visible effects on the position of 
the doppler peaks, which are due to the oscillatory phase 
with which the photon-baryon fluid meets the last scatter- 
ing band (see Hu & Sugiyama 1995). The photon-baryon 
fluid oscillates as (xis{krs), where k is the comoving scale 

and rs is the sound horizon (rs = /J^^"^ dTCs{a), Cs(a) is the 
sound speed). Given the photon-baryon ratio, rs{a) follow a 
similar trend as r(a). Since in volatile models r{a) is smaller 
than in CDM, so will be rs{a), and the peaks of the spec- 
trum will appear in correspondence to higher k ( i.e. higher 
I) values. 

This is a specific features of these models, in neutrino 
models the same effect plays a role, but shifting the peaks 
in the opposite direction (Dodelson et al. 1996). 

For a given n, the height of the peaks is fixed by (i) 
the ratio between baryon and photon densities, i.e. Q^h^, 
and (ii) the ratio between matter and radiation densities. 
At fixed Qi, and h the main reason for a higher doppler peak 
in volatile models (with respect to CDM) is the delayed 
matter-radiation equivalence, for which both SMLC and, 
possibly, volatilcs can be responsible. In neutrino models 
without ad-hoc SMLC, only the possible delay due to late 
derelativizing v's may exist. This is why volatile Ci spec- 
tra and standard neutrino ones look so different. However, 



there is a tiny further contribution in the boost of the peak 
due to the free-streaming of the hot component. Several au- 
thors (Ma & Bertschinger 1995, Dodelson et al. 1996) have 
shown that even in high Zder neutrino models the doppler 
peaks arc enhanced with respect to CDM, and in that case 
the free-streaming of the hot component is to be considered 
responsible for the enhancement. Free-streaming, in fact, 
causes a decay in the potential $ which contributes as a 
forcing factor (trough "l?) in the equations whose solution are 
the sonic oscillations in the photon-baryon fiuid, displacing 
their zero-point and, henceforth, the phase by which they 
enter the last scattering band. In the standard neutrino case, 
this effect causes a variation of 10% at most on the C;, and 
typically of 2% on the first doppler peak. 

In principle one can expect that the different momen- 
tum distribution of volatiles may alterate the free-streaming 
behaviour. Such differences, if they exist, can be found by 
comparing volatile spectra with the technical neutrino ones. 
The differences between the two spectra are presented are 
shown in fig. 4, and amount to 2% at most. Although mod- 
est, this is another feature that characterizes volatile models 
with respect to neutrino one. 

In comparison with such finely tuned predictions from 
theoretical models, currently available data are still affected 
by huge errorbars. However, some feature seems already 
evident from them. In fig. 6-15 we perform a comparison 
of model predictions with data and show that the doppler 
peak observed by the Saskatoon experiment (Netterfield et 
al. 1997) exceeds the one expected in pure CDM once it is 
normalized to COBE data (Bennet et al. 1996). While it 
is evident that volatile models show a higher doppler peak. 
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Figure 4. differences in the radiation power spectrum between 
the volatile case and the techical neutrino case. The model pa- 
rameters are the following: 1) f2^ = 0.16, d = 4.25; 2) = 
0.22, d = 1.03; 5;f2h = 0.22, d = 2.06; .^j fi^ = 0.33, d = 2.06. 



it is clear that a fit could be reached also changing other 

parameters, e.g., by taking n > 1. In fig. 4 we show what 
happens in neutrino models if the spectrum is anti-tilted to 
n — 1.1 and to n = 1.2. Indeed, the first doppler peak is 
raised (which is desirable), but also the following peaks are 
raised, making difficult the agreement with the results from 
the CAT experiment (Scott et al. 1996). In section 5 similar 
considerations will be used in order to constrain the whole 
set of volatile models. 



4 CONSTRAINTS FROM LARGE SCALE 
STRUCTURE 

Mixed model parameters can be constrained from particle 
physics and/or from LSS. In this section we review a num- 
ber of the latter constraints, which can be tested without 
discussing non linear evolution. In Appendix A we debate 
constraints on the spectral index n arising from inflation. 

Even without considering their non-linear evolution, 
models can be constrained through the following prescrip- 
tions: 

(i) The numerical constant A<s, in the spectrum (1.1), 
must give a value of Ci consistent with the COBE 
quadrupole QrmB,ps. Values of consistent with the 
Qrms,PS values, for a given n, within 3cr's, can be kept. 

(ii) COBE quadrupole therefore fixes the normalization 
at small k. The first large k test to consider, then, is the 
behaviour on the 8 h~^Mpc scale. Quite in general, the mass 



o 




Figure 5. a standard neutrino model with different n values. Ci 
curves from bottom to top correspond to n = 1, 1.1, 1.2, 1, 3. Bold 
data points refer to COBE (Tegmark 1996), CAT (Scott et al 
1996), Saskatoon (Netterfield et al. 1997); see Scott et al (1995) 
for a summary of all the experiments. 

M\, within a sphere of radius L — A/i^^Mpc, is 

Mx = 5.96 ■ lQ'^Q.h''MQ{\/?,f . (4.1) 

Therefore, the 8/i~^Mpc scale is a typical cluster scale. 
Here optical and X-ray data are to be exploited to work 
out the mass variance as and models should fit such ob- 
servational outputs. Optical data provide the cluster mass 
function through a virial analysis of galaxy velocities within 
clusters. X-ray determinations, instead, are based on obser- 
vational temperature functions. If clusters are substantially 
virialized and the intracluster gas is isothermal, the mass M 
of a cluster can then be obtained, once the ratio f3th,gai be- 
tween thermal or galaxy kinetic energy (per unit mass) and 
gravitational potential energy (per unit mass) is known. Val- 
ues for /3's are currently obtained from numerical models. 

Henry & Arnaud (1991) compiled a complete X-ray 
fiux-limited sample of 25 clusters which is still in use for 
such determinations. Assuming an isothermal gas, full viri- 
alization and f5 = figai/Pth = 1-2 they had estimated 
Us = 0.59 ± 0.02. Their error does not include /3 uncertainty. 
Various authors then followed analogous patterns (see, e.g.. 
White et al. 1993, Viana & Liddle 1996). Recently Eke et 
al. (1996) used Navarro et al. (1995) cluster simulations to 
take /9 = 1 with an error ;$ 6 %. Accordingly they found 
as = 0.50 ± 0.04. By comparing the above results one can 
estimate that, to obtain as ~ 0.7, under the assumption 
of full virialization and purely isothermal gas, j3 ~ 1.4 is 
needed. 

An estimate of cluster masses independent from clus- 
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ter models can be obtained by comparing optical and X-ray 
data. Recent analyses (Girardi et al. 1998) seem to indicate 
values of /3 ~ 0.88. In our opinion, such outputs do not 
strengthen the case of a safe cluster mass determination, 
as they arc more than 12 % below Navarro ct al. (1995) 
ratio and might indicate a non equilibrium situation. Fur- 
thermore, it ought to be outlined that cluster mass determi- 
nations based on a pure virial equilibrium assumption con- 
flict with the observed baryon abundances and would require 
cosmological models with Q.b ~ 0.16-0.20, in contrast with 
BBNS constraints, if all dark matter is CDM and f2 = 1. 

If HDM is only partially bound in clusters and their 
masses arc underestimated by ^ 15-20 %, the latter conflict 
can be overcome. (Alternative way outs, of course, are that 
n < 1 or A 7^ 0.) Therefore, in order that data be consistent 
with mixed models, some mechanism should cause a slight 
but systematic underestimate of cluster masses. Owing to 
such uncertainties, wo can state that cluster data constrain 
as, within the interval 0.46-0.70. 

These constraints can also be expressed with direct 
reference to the cumulative cluster number density. Defin- 
ing the meiss M{u) for which the top-hat mass variance 
f"M(u) = 5c/u (here 5c values from 1.55 to 1.69 can be con- 
sidered) the Press & Schechter approach yields the number 
density 

/>oo 

n(> M) = ^j2l^{plM) \ dw[M/M(u)] exp(-uV2) . 

(4.2) 

A usual way to compare it with data amounts to taking 
M = 4.2ft"^ ■ IO^'^Mq and considering then Nd = n(> 
M)(100/i"^Mpc)^ for the above M value. With a range of 
uncertainty comparable with the one discussed for cts, op- 
tical and X-ray data converge towards a value of N^i ~ 4. 
Henceforth viable models should have 1^ N^i^ 10, for one 
of the above values of <5c . 

There is a slight difference between testing a model in 
respect to crs or N^i- This amounts to the different impact 
that the slope of the transferred spectrum has on expected 
values. Observations, however, also constrain the observed 
spectral slope, as we shall detail at the point (iv). 

(iii) In order to have A^ci and a% consistent with obser- 
vations, the interval obtained from COBE quadrupole 
may have to be restricted. The residual range of A-^ values 
can then be used to evaluate the expected density of high 
z objects, that mixed models risk to under-produce. The 
most restrictive constraints comes from computing figas = 
aQ.h^co\\ in damped Lyman a systems (for a review see 
Wolfe, 1993). 

It can be shown that 

acoii = erfc[<5c/V2cr(M, z)\, (4.3) 

where cr(M, z) is the (top hat) mass variance (for mass M 
at redshift z) and a is an efficiency parameter which should 
be i$ 1. More specifically, using such expression, one can 
evaluate DLas{M, z) = figas x 10^/a. Then, taking z = 4.25, 
5c = 1.69 and M = 5-10^ h~^MQ we have a figure to compare 
with the observational value given by Storrie-Lombardi et 
al. (1995): Dlus = 2.2 ± 0.6. 

Only models for which the predicted value of Dlus ex- 
ceeds 0.5, at least for a part of the allowed A<sf interval, are 



therefore viable. In turn, also for viable models, this may 
yield a further restriction on the interval. 

(iv) Models viable in respect to previous criteria should 
also have a fair slope of the transferred spectrum. Its slope 
can be quantified through the extra power parameter F = 
7.13-10"^ (cr8/cr25)^°'^'^ ((^8,25 are mass Variances ou the scalcs 
R = 8,25/i-^Mpc). Using APM and Abell/ACO samples 
Peacock and Dodds (1994) and Borgani et al. (1997) ob- 
tained F in the intervals 0.19-0.27 and 0.18-0.25, respec- 
tively. Such intervals essentially correspond to 2 a's. Further- 
more the lower limit can be particularly sensitive to under- 
estimates of non-linear effects. Henceforth, models yielding 
F outside an interval 0.13 0.27 are hardly viable. 

One can also test models against bulk velocities recon- 
structcxi POTENT from observational data. This causes no 
constraint, at the 2 a level, on models which survived previ- 
ous tests. 

In BP a number of plots of the transferred spectra of vi- 
able models, were shown against LCRS reconstructed spec- 
tral points (Liu ct al. 1996). However, previous constraints 
include most quantitative limitations and models passing 
them fit spectral data. Fig. 1 is taken from BP and reports 
the curves on the fl, d plane limiting areas where viable 
mixed models exist for various primeval n values, if fii, = 0.1. 
All models considered in the next sections, both for Oj, =0.1 
and fib = 0.05, were previously found to satisfy the above 
constraints. 



5 CMB DATA AND PARAMETER SPACE 
LIMITATION 

In this section we give the CBR spectra of the hot-volatile 
models and compare them with available data, ranging from 
l = 2tol< 500. 

We evaluated the spectra for several parameter choices 
allowed by LSS constraints (see fig. 1). Significant example 
of Cj spectra are shown in figs. 6-15 while the corresponding 
LSS predictions axe summarized in table 1. 

Some parameter sets are compatible with neutrino hot 
dark matter, while models with a low Zder and n> 1 arc ob- 
tainable with volatile hot dark matter. Since the height of 
the first doppler peak is very sensitive to the baryon abun- 
dance, we considered two values of Q.b, namely 0.05 and 0.1. 
Spectra are normalized to Qrms,ps assuming no contribu- 
tion of gravitational waves. As is known, their contribution 
would raise the low-1 tail of the C; spectrum, therefore re- 
ducing the gap between the Sacks- Wolfe plateau and the top 
of the first doppler peak. 

Models with Qb = 0.1 systematically show a peak less 
pronounced than models with fij, = 0.05. It is well known 
that models with a given h and hot component show a lower 
doppler peak for smaller f2(,; in top of that, here there is 
a further effect: LSS constraints often arc compatible with 
a part of the observational Qrms,ps interval, and low Q,b 
models tend to be consistent with low Qrma,ps values. 

While in fig. 5 we plot most data are available, in figs. 6- 
15 we compare models with data from COBE (Tegmark 
1996) Saskatoon (Netterflled et al. 1997) and CAT (Scott 
et al. 1996) only. 

Figs. 6-15 show a systematic trend: for a given large / 
value, Ci increases with both n and d = IQ"^ /zder- On the 
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[h] 

Table 1. Model parameters and power spectra. Column 1: 
volatile fractional density; Column 2: redshift at which the volatile 

component becomes non-relativistic {d = lO^/z^er); Column 3: 
total number of C(]nivalcnt masslcss neutrinos (n^ = g^,^j j 
Column 4: n value considered; Columns 5-7: Large scale structure 
predictions. 
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Figure 6. volatile spectra are compared with the observational 
data from COBE, Saskatoon and CAT experiments. The error- 
bars correspond to 1 a errors. Solid bold lines refer to standard 
CDM (h = 0.5, n = 1, Q = 1) and a neutrino model with Qh = 
0.3. Shaded areas are the volatile spectra with parameters speci- 
fied in the figure. Upward shading refers to f2 = 0.1 models while 
downward shading to Q = 0.05 Notice that z^i^j. = lO^'/d = 1250. 
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contrary, for a given large k value, the matter fluctuation 
spectrum P{k) increases with n but is damped for large d, 
so that these to effects tends to compensate. 

This is one of the reasons why LSS constraints can be 
compatible with n as high as 1.4. On the contrary, figs. 14- 
15 show that CBR spectra already disfavour n = 1.3 if 
d> 4 («der^ 2500) is considered, no matter the value of Q.b. 
Volatile models with n > 1.2 are largely out of the errorbars, 
and should be considered as scarcely viable. Nevertheless, 
even for n = 1.2, volatile models allow a higher first doppler 
peak without raising the following ones, and therefore fit the 
data better than neutrino models. Just as large n, also large 
d causes conflict with data, by itself. For example, fig. 10 
show that models with d = 16 are disfavoured, even with 
low fih and n — 1.1. 

As pointed out in section 2, volatile models require a 
sterile component whose energy density is proportional to 
Qh d. Its effective number of degrees of freedom is linked 
to the equivalence redshift, which in turn affects both the 
shape parameter F and the height of the first doppler peak. 
Dodelson et al. (1994) considered the matter power spec- 
trum in the case of a r-neutrino decay (rCDM model), and 
found that even in that case the effective number of degrees 
of freedom (/,y,e// is bigger that in standard CDM. They out- 
lined that, in order to lower F at least down to 0.3 (Peacock 
& Dodds 1994), in a /i = 0.5 universe with n = 1, an equiv- 
alent number of massless neutrinos as high as 16 is needed. 



White et al. (1995), who also consider tCDM models but 
with a lighter neutrino, also pointed out that the predicted 
F of these models is lower due to the high gv,eff, and show 
that a lower F implies a higher first doppler peak. Their 
work, however, is only qualitative, and they don't infer any 
restriction in the parameter space using the data. Lookng 
at the data, we found out that if n = 1.1, CBR data mod- 
els with an equivalent number of neutrino species N^> 10 
as in figs. 10, 12 and 15 are disfavoured . Models like the 
one shown in fig. 8 [Ni, ~ 7) seem to better fit the data, 
although even lower N^, (~ 5), as provided by the model in 
fig. 9, should not be disreguarded. 

Keeping to n = 1, LSS already exclude very high Qh d 
values, so that a low is automatically ensured. The mod- 
els shown in figs. 6-7 seem to well fit the data, with a cor- 
responding Nv = 5 — 7. 



6 CONCLUSIONS 

In this work we have analizcd mixed models from the point 
of view of both LSS and CMB predictions. We considered 
different hot dark matter components: the standard neutrino 
case and the volatile case in which particles come from the 
decay of heavier ones. 

First we tested the mixed models on available LSS data 
requiring fair predictions for as, F, DLAS and Nd. This 
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Figure 15. volatile spectra compared to the 
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analysis shows that it must be 0.3. This comes as no 
surprise, as mixed models with greater have not been 
considered since long. The new result is that taking n up 
to 1.4 does not ease the problems previously found for large 

On the contrary, volatile models together with n > 1 
significantly widen the parameter space in the low Zder direc- 
tion and viable models even with Zder ~ 600 can be found. In 
fact, as far as P{k) is concerned, we found a nearly degener- 
ate behaviour of the parameters n and Zder, as the damping 
on the high k values due to low Zder can be compensated by 
high n. 

CBR data, apparently, break the degeneracy. In sec- 
tion 3 we have shown that the CBR spectrum of volatile 
models is significantly different from standard CDM and also 
from neutrino models usually considered. In fact SMLC and 
late Zder volatiles cause a late z^q and, henceforth, a higher 
first doppler peak. Minor effects are caused by the typical 
momentum distribution of volatiles. These effects amounts 
to 2 % at most in the C; spectrum and only an accurate anal- 
ysis of the results of future satellites, as MAP and Planck, 
could allow to detect it. 

CBR spectra of volatile models were then compared 
with available data from different experiments, namely those 
from COBE, Saskatoon and CAT experiments. Data on 
CBR spectrum at large I imply that temperature fluctua- 
tions AT/T ~ 10~^// are appreciated. Therefore, measures 
of the CBR spectrum, for high / values, still need to be 
treated with some reserve. It seems however clear that re- 
cent observations tend to indicate a doppler peak higher 
than expected both for pure CDM and for mixed models 
with early derelativization, such as most neutrino models. 



Taking n > 1 and/or late derelativization raises the doppler 
peak and affects the CBR spectrum at high I. The first ques- 
tion we tried to answer is how far we can and have to go from 
pure CDM and n = 1 to meet current large / data. 

We found that volatile models could cure this discrep- 
ancy, while ensuring a viable scenario for structure forma- 
tion. 

In turn, large / data imply restrictions in the parameter 
space, complementary to the ones derived from LSS while a 
fit of such data requires only a slight departure from pure 
CDM and n — 1. allows us to say that mixed models are 
in very good shape. For example, fig.s 8-9 show the Ci be- 
haviour for n = 1.1 and HDM ranging from 11% to 16 %. 
Such models provide excellent fits to current data and, as 
explained in BP, are also in agreement with LSS. 

Other models, for larger n and Qh or lower Zder, show 
only a marginal fit with current observations. Hopefully, fu- 
ture data on high Vs will be more restrictive and allow safer 
constraints. At present, such models cannot be ruled out, al- 
though they are more discrepant from pure CDM and n = 1 
than high I data require. 

In our opinion, however, CBR data can already be said 
to exclude a number of models which fitted LSS data. In 
general, models with n > 1.3 and Zder < 1000 seem out of 
the range of reasonable expectations. 

Altogether, three kinds of departures from CDM and 
Zel'dovich were considered in this work: large Qh, low Zder 
and n > 1. Large (but allowed) ilh values, by themselves, do 
not ease the agreement of models with high I data. Taking 
n > 1 eases the agreement of models with data for I ~ 
200, as is expected, but seems to rise the angular spectrum 
above data for greater Vs. Taking low Zder, instead, raises the 
doppler peak, but does not spoil the agreement with greater I 
data. Current data, therefore, seems to support models with 
a limited amount of HDM or volatile materials, possibly in 
association with n slightly above unity, to compensate some 
effects on LSS. 

Note that the analysis of this work is carried out keep- 
ing h — 0.5, allowing for no cosmological constant, and con- 
straining the total density to be critical. E.g., raising h would 
probably allow and require a stronger deviation from pure 
CDM and n = 1. We plan to widen our analysis of the pa- 
rameter space in the near future, also in connection with 
the expected arrival of fresh observational data on the CBR 
spectrum. 
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7 INFLATIONARY MODELS YIELDING N > 1 

This section is a quick review of results in the literature, 
aiming to show that there is a wide class of inflationary 
models which predict n > 1, but ^ 1.4. 

During inflation, quantum fluctuations of the inHaton 
field (fi on the event horizon give rise to density fluctuations. 
Their amplitude and power spectrum are related to the Hub- 
ble parameter H during inflation and to the speed (f of the 
slow-roUing-down process. The critical quantity is the ratio 
W{k) = H'^/ifi, where H and ip axe taken when the scale 
2n/k is the event horizon. It can be shown that 



1 + 2 



d(log W) 
d(log k) 



(a.l) 



and, if W (slowly) decreases with time, we have the stan- 
dard case of n (slightly) below unity. Such decrease occurs if 
the downhill motion of ifi is accelerated and an opposite be- 
haviour occurs if ifi decreases while approaching a minimum. 
The basic reason why a potential yielding such a behaviour 
seems unappealing, is that the very last stages of inflation 
should rather see a significant i/p-field acceleration, ending 
up into a regime of damped oscillations around the true vac- 
uum, when reheating occurs. 

However, the usual perspective can be reversed if the re- 
heating does not arise when an initially smooth acceleration 
finally grows faster and faster, but is triggered by an abrupt 
first order phase transition, perhaps due to the break of the 
GUT symmetry. Before it and since the Planck time, most 
energy resided in potential terms, so granting a vacuum- 
dominated expansion. This picture of the early cosmic ex- 
pansion is the so-called hybrid inflation, initially proposed 
by Linde (1991a). 

A toy-model to realize such scenario (Linde 1991b, 
1994) is obtainable from the potential 



(a.2) 



which depends on the scalar fields ip and x, that we shall 
see to evolve slowly and fastly, respectively. If the slow field 
is embedded in mass terms, the potential reads 



(a.3) 



where 



a4 4 , 2 2 

A — fj, + m tp 



and 



= 2{g'^</ - \^?) . (a.4) 



Eq. (a.3) shows that V has a minimum at x = 0, provided 
th at > 0. If < 0, instead, the minimum is for x~ 
y^MV2A, yielding ^ when = 0. 

Large ip values therefore require that x vanishes and 
then the potential 



V{p>, 0) = + m^<p^ 



(a.5) 



gives a Planck-time inflation, as V goes from an initial value 
~ tpf to a value ~ /i*. The downhill motion of tp will de- 
celerate as soon as the second term at the r.h.s. of eq. (a.5) 
becomes negligible, in respect to which acts as a cosmo- 
logical constant. This regime breaks down when the critical 
value ipc ~ \/A/i/(? is attained. Then changes sign and 
the configuration x = becomes unstable. We have then a 
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transition to the true vacuum configuration x, which reheats 
(or heats) the Universe. 

There are several constraints to the above picture, in 
order that at least 60 e-foldings occur with (p > and 
fluctuations have a fair amplitude. They arc discussed in 
several papers (see, e.g., Copeland et al. 1994, and references 
therein) and cause the restriction 1.4. 

It is fair to outline that hybrid inflation is not just one 
of the many possible variations on the inflationary theme. 
In spite of the apparent complication of the above scheme, it 
is an intrinsically simple picture and one of the few patterns 
which can allow to recover a joint particle-astrophysical pic- 
ture of the very early Universe. 
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